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Correspondence between classical periodic orbits and quantum shell struc- 
ture is investigated for a reflection-asymmetric deformed oscillator model as 
a function of quadrupole and octupole deformation parameters. Periodic or- 
bit theory reveals several aspects of quantum level structure for this non- 
integrable system. Good classical-quantum correspondence is obtained in the 
Fourier transform of the quantum level density, and importance of periodic 
orbit bifurcation is demonstrated. Systematic survey of the local minima of 
shell energies in the two-dimensional deformation parameter space shows that 
prominent shell structures do emerge at finite values of the octupole param- 
eter. Correspondences between the regions exhibiting strong shell effects and 
the classical bifurcation lines are investigated, and significance of these bifur- 
cations is indicated. 
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1. INTRODUCTION 



Shell structure is one of the important aspects of finite quantum many-body systems. 
In the single-particle level density, one may generally find some regular patterns like shells 
consisting of dense and thin regions. This pattern changes with deformation, and the system 
favors the shape which makes the level density at the Fermi surface lower. Predictions of 
the superdeformed (extremely large quadrupole deformations with axis ratio about 2:1) and 
the hyperdeformed (the axis ratio about 3:1) nuclear states, which are hot current topics 
of high-spin nuclear structure physics, had been based on the above kind of consideration. 
Strong shell effects which appear at ellipsoidal shapes with the axis ratios about 2:1 and 
3:1 play essential roles in stabilizing such exotic shapes. Reflection-asymmetric degrees 
of freedom are also one of the most exciting subjects in the current high-spin physics. 
Superdeformed potentials possess remarkable single-particle level structures where levels 
with different parities approximately (exactly in the harmonic oscillator limit) degenerate 
in the same major shell, which may bring about strong octupole correlations ||l|-^. Recent 
remarkable development of large 7-ray detector arrays encourages us to find such exotic 
shapes like reflection-asymmetric superdeformations, and to investigate the mechanism of 
producing them. Recently micro-cluster physics have also attracted much attentions, and 
many nuclear physicists have been contributing to this new field. Shell structures and 
deformations of clusters are very interesting subjects — their shapes can actually be seen 
with an electron microscope — , and one can apply almost the same theoretical framework 
to both nuclei and micro-clusters 

A clear understanding of the origin of shell structure may be obtained by the use of 
semiclassical theory. Correspondences between quantum spectra and classical dynamical 
properties of Hamiltonian systems have been extensively investigated for two limiting cases, 
namely, for integrable and strongly chaotic situations. Most physical systems are situated 
in the midst of these limits, however, and belong to what we call 'mixed' systems. The 
semiclassical theory for mixed systems is difficult and only few aspects have been clarified 
up to now. This difficulty is associated with the periodic-orbit bifurcations (characteristic to 
the mixed systems) where the stationary phase approximation (SPA) and the conventional 
trace formula for representing the quantum spectrum in terms of classical periodic orbits 
breaks down. Fortunately, however, an approach in the inverse direction sometimes works 
and one can extract the periodic orbit information from the quantum spectrum by means of 
the Fourier transformations. This approach is very useful to understand the shell structure 
of the quantum spectrum. We take this approach and clarify some aspects of a mixed 
system, directing our attention to the influence of the bifurcations of short periodic orbits 
on the gross structure of the quantum spectrum. 
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In this paper, we investigate the classical-quantum correspondence for an axially- 
symmetric deformed oscillator model with reflection-asymmetric terms. This is a non- 
integrable model and chaotic behavior gradually emerges in the dynamics as the octupole 
deformation becomes large. The role of periodic orbit bifurcations will be emphasized. In 
section |[ basic elements of the semiclassical theory relevant to our analysis are briefly re- 
viewed. Special attention will be paid on the classical bifurcation phenomena and their 
effects on the quantum spectra. In section ^ our model is introduced and several aspects of 
it are summarized. In sections ^'^-j^, we will present numerical results of the semiclassical 
analysis and discuss their implications. It will be shown that prominent shell structures 
emerge for finite octupole deformations superposed on the prolate shapes. Origins of such 
new shell structures will be clarified. Section H is devoted to summary and conclusion. 



2. BASIC FORMULAE 

2.1. Classical Hamiltonian system 

Let us consider a Hamiltonian system with / degrees of freedom. The equation of motion 
(EOM) for a 2/-dimensional phase space vector Z = {p,q) is expressed as 

4-Z = A^, ^=1 ; ), (2.1) 




dt dZ' 

where O and / denote the /-dimensional zero and identical matrices, respectively. Now 
consider a bundle of trajectories around a certain solution Za{t) and write them as Z{t) — 
Za{t) + 5Z{t). Then the EOM for 6Z{t) is given by 

±8Z^AniZ^it))5Z, H(Z).,=g| (2.2) 

up to the first order in 6Z . H is called Hessian matrix. One can easily integrate the above 
differential equation and obtain the following solution: 



SZit) = exp 



[ drAHir) 



dZ{to) = Sa{t-to)SZ{to), (2.3) 



where Tr denotes that the exponential is defined by time-ordered product. Sa is called 
stability matrix of the trajectory a, whose eigenvalues determine its stability. 

Let S denote a (2/ — 2) dimensional hypersurfacc in the phase space with fixed energy 
E. It defines a time-discretized mapping M : S S with classical trajectories, which is 
called the Poincare map. Periodic orbits Z are defined as the fixed points of A4, namely, 
by M{Z) = Z. The linear part Mr (with respect to SZ) of A4 about a periodic orbit Zr is 
called "monodromy matrix" and describes the stability of the orbit: 

M{Zr + 6Z) = Zr + Mr6Z + Oi6Z^). (2.4) 

The monodromy matrix is a symplectic matrix satisfying 

AMJ A-'^ = Mr^^, (2.5) 
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and this property restricts its eigenvalues as follows. Let A be one of the eigenvalues of Mr- 



Relation (2.5) guarantees that the reciprocal of A is another eigenvalue of Mr- Furthermore, 
Mr is a real matrix so that the complex conjugates of these eigenvalues are also eigenvalues. 
Thus the eigenvalues of the monodromy matrix generally appear in quartets (e^"^'*^). In 
the two-dimensional case, they appear in a pair. 

Let us now proceed to a discussion on bifurcations of stable periodic orbits. Consider a 
trajectory that emerges at q with energy E and returns to the initial point, q' = q. This 
kind of trajectory certainly exists for any q. The condition for this trajectory to be periodic 
is that the initial and final momenta coincide; namely, p' = p. Using the Hamilton- Jacobi 
equation, this condition can be rewritten as 

O^P'-P^ - (- = M^, (2.6) 

oq' \ oq J oq 

where S{q) denotes the action integral along the closed path under consideration. Thus, 
the periodic orbit is the stationary point of this action integral. Let us expand S{q) about 
the periodic point q: 

= S+^6q^B5q + ---. (2.7) 

After some simple matrix rearrangements, one can express B in terms of the quadrants of 
the monodromy matrix as 



B = B-{I-A)C-^(I-D), M=\ \. (2 




Note that 

det(l -M) = - dct[C] det[S]. (2.9) 

The above relations provide us with a clear understandings of the connection between eigen- 
values of the monodromy matrix and the bifurcation of periodic orbits. Suppose that one 
of the eigenvalues of M becomes unity. Then the curvature tensor B for the action 5" has a 
zero eigenvalue. This means that the stationary points of S (periodic orbits) locally form a 
continuous set and a bifurcation can occur hereafter; namely a (few) new stationary point(s) 
can emerge. 



2.2. Trace formula 

By means of the semiclassical theory, we can relate properties of the quantum spectrum 
with those of the corresponding classical system. For non-integrable Hamiltonian systems, 
the Gutzwiller's trace formula represents the quantum level density g{E) = 5{E — 
Ei) as a sum over classical periodic orbits: 

g(^^\E)^-g{E) + Y.^r.r{E)cos _ ^\ ^^.iq) 
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g{E) is called Weyl term (or Thomas- Fermi approximation), which is a monotonic function 
of energy. The sum is taken over all primitive periodic orbits r and their multiple traversals. 
Sr — §j.P ' dq is the action integral along the orbit r, and fir denotes the Maslov phase. 
Several numerical application of this formula to strongly chaotic systems have shown its 
effectiveness. As is well known, however, exact reproduction of a quantum spectrum is not 
an easy task, because one has to treat huge number of periodic orbits which show exponential 
proliferation as a function of energy. On the other hand, since our purpose is to understand 
the gross structure of the spectrum we need only finite number of periodic orbits which have 
rather short periods. If one is interested in the gross structure with energy resolution 5E^ 
the change of the phase in Eq. (2.1C) must be less than 27r for the change of energy by 5E. 
Namely, 

nSr{E + 5E) - nSr{E) ~ n !r ' 5E = nTrSE < 27rh, (2.11) 

oE 

which leads to the relation 

nTr<T^e..= ^. (2.12) 

Thus, we need only short periodic orbits whose periods are less than Tmax defined above. Al- 
though much efforts have been devoted to reproduce individual eigenenergies by calculating 
millions of periodic orbits, gross structures of the level spectra have rarely been discussed in 
connection with periodic orbits for non-intcgrable systems. (For integrable systems, there 
are several works; see for instance, Refs. flp^-p^ ) 



2.3. Bifurcations 

Let us next discuss the condition for the amplitude factor Anr in the trace formula ( |2.10| ) to 
take a large value. In the stationary phase approximation, the amplitude factor is expressed 
for isolated orbits as 

1 Tr 

^ ^v/|det(l-lw7)I' ^^'^^^ 



For degenerate orbits with degeneracy 1, it is represented by 



An Br '■^'^ 



i2nn)y^ ^|det(l -M/)r 



Br= dt 







dipjt + Tr) 
dp^t) 



-1/2 



(2.14) 



where (p and p^p denote an ignorable variable in the Hamiltonian and its canonically conjugate 
momentum, respectively (see Appendix in Ref. [Q). From the difference in order of h, one 
sees that the contribution of degenerate orbits is more important than that of isolated ones. 

Another important factor, which plays an essential role in our analysis below, is the stabil- 
ity factor det(l — Af,."). Its value is independent of the point chosen on the periodic orbit. As 
discussed above, eigenvalues of the monodromy matrix Mr appear in pairs (+/— )(e'^, e^"^), 
A being real or pure imaginary, or in quartets (e*"^*''). One should note that the periodic 
orbit generally appears in at least one parameter family, so that the monodromy matrix 
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always has two unit eigenvalues. Other pairs of unit eigenvalues correspond to the global 
continuous symmetries which the Hamiltonian possesses but the orbit itself does not. These 
degrees of freedom are responsible for the degeneracy and, as seen in Eq. ( ^.14 ), can be 
separated out from the definition of Mr appearing in the stability factor. In the case of 
three-dimensional systems with axial symmetry, the monodromy matrix has (except for iso- 
lated non-degenerate orbits) four unit eigenvalues, and the remaining two appear in a pair 
(-1-/— )(e'^, e~^). Thus, may be reduced to a (2x2) matrix for most orbits. A is purely 
imaginary for stable orbits, and real for unstable ones. The stability factor for each case 
becomes 

det(l - Mr) = 2 - TiMr 

with 

TrMr^l^''''^^\ ^ = (2.15) 

I ±2cosh(a)) unstable, \ — a 

Thus, the stability of orbit is determined by the value of TrM. 

If a parameter in the Hamiltonian is continuously varied, the periodic orbits change their 
shapes and the values of A also change continuously. It may occur that the f3 for a certain 
stable orbit becomes a fraction of 27r, namely, f3 = 27rTO/n with n and m being relatively 
prime integers. At this point TrM^" becomes 2, and the amplitude factor Anr suffers 
divergence. This singularity corresponds to the period n-upling bifurcation of the orbit r. 
Near the bifurcation point, the stationary phase approximation breaks down. It is then 
necessary to take into account higher-order fluctuations about the classical orbit to extract 
a finite value of Anr fl^]. Although we leave this task as an challenging future subject, we 
expect that the amplitude factor takes large value in the bifurcation region. It will result 
in a large-amplitude oscillation in the level density, leading to an enhancement of the shell 
effect. 



3. THE MODEL AND ITS SCALING PROPERTIES 
3.1. The model 

We adopt a model Hamiltonian consisting of an axially deformed harmonic oscillator and 
a reflection asymmetric octupole deformed potential: 

= ^ + E - ^3oM-o^ [r'Y,om " . (3.1) 

i 

Here, the double primes denote that the variables in square brackets are defined in terms 
of the doubly-stretched coordinates x" = {uJi/u;o)xi, where = i^^xi^yi^zY^^ being deter- 
mined so that the volume conservation condition is satisfied fl^ . For simplicity, we define 
dimensionless variables as follows: 



Pi — > ^/Mhojo Pi 

Qi — > y^h/Mujaqi (3.2) 
H — > hcjoH 
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Then, the Hamiltonian is written as 
H = -p 



(3.3) 



2" 2 

Since the radial dependence of the octupole potential is quadratic, the Hamiltonian obeys 
the following scaling rule: 



H{ap,aq) = a^H{p,q). 



(3.4) 



Thanks to this property, if one solves the Hamilton equations of motion at a certain energy 
i?o, the solution at any energy E is obtained by just scale transforming the solution at Eq 
as Z{t\ E) — y/E/Eo Z{t; Eq), Z denoting a phase space vector (p, q). 



In the cylindrical coordinates (p, z, ip), the above Hamiltonian (3.3) is written as 



P^ 
2p2 



+ 2^ / 7 2z^ - 3zp^ 

-z A30 \ 



(3.5) 



Thus, we can treat the system as a two-dimensional one with fixed angular momentum p^. 
This reduction enables us to make use of the Poincare surface of section in order to survey 
the classical phase space profile. Note that x = P^pl E is a scaling-invariant parameter and, 
therefore, classical properties are the same for the same x- 

In the following sections, we shall investigate how the properties of the quantum spectrum 
for the Hamiltonian (3.1) changes as the two deformation parameters, 5osc = (^± — ^z)/^ 
and A30, are varied. We shall then discuss the physical origins of these changes by means of 
the periodic orbit theory reviewed in section ^. 



3.2. Fourier transformation of quantum level density 

As we will see in the following numerical analyses, the Hamiltonian above becomes chaotic 
with increasing octupole deformation parameter A30. But considerable parts of the phase 
space remain regular and the system is considered as a so-called 'mixed system'. The trace 
formula based on the stationary phase approximation (SPA) does not work well in such 
situations. The amplitude factors suffer divergences at the bifurcation points of stable 
periodic orbits because of the breakdown of SPA. Consequently, we cannot directly use the 
conventional semiclassical expression to analyze the quantum spectrum. Fortunately, we can 
avoid the above difficulty by using the Fourier transformation technique for the quantum 
level density. Suppose that the level density is characterized by the classical periodic orbits 
and is expressed as 

9^i^ -EE Anr{E) COS - |/.„.) , (3.6) 

n— 1 r ^ ^ 
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without specifying the concrete expression of the amphtude factor A„r{E) which may be 



obtained by going beyond the SPA. Since our model obeys the scahng rule (3.4), energy 



(3.7) 



dependence of the classical variables entering Eq. (3_^) is factored out as follows: 
Sr{E) — E ■ Tj., 

A„,(£;) = i?'^'V2Ai°,\ 

where dr is the effective degeneracy of the orbit, and equal to 1 for most orbits due to the 
axial symmetry. For isolated orbits, = and we expect that their contributions to the 
level density may be small. The degeneracies are integers in the classical dynamics, but 
in the quantum mechanics this restriction is relaxed and dr changes continuously in the 
bifurcation regions Using the above relations, Eq. ( ^77| ) is expressed as 

)(i?)=5(i?)+^i?'^'-/^4°)cosf-"^^^ 



V ^ 



2 Mnr 



(3.8) 



Now let us consider the Fourier transformation of the level density, defined by 

2' 



F(s) dE e'^'^'^ E-'^/^g(E) ■ exp 

27rn, J 



E 



E„ 



(3.9) 



Here the Gaussian damping factor is used for energy cut-off, and we shall put d = 1 in order 
to cancel the energy dependence of the amplitude factors for most orbits. Inserting the 
quantum level density g{E) = 5{E — Ei) and the semiclassical one (3.8) into Eq. (^j|), 
we obtain quantum mechanical and semiclassical expressions for F{s): 



/E, 



F(-)(.)c.F(5) + ^A(°) 



exp 



E, 



1 / s — nTr 
2 



As 



(3.10) 



(3.11) 



respectively, where As = h/E^i^^. F^^^™^ is calculated from the single-particle spectrum 
obtained by diagonalizing the Hamiltonian with deformed oscillator basis. The result is 
compared with the semiclassical expression F^^'^K In Eq. (3.11), F{s) corresponds to the 
Weyl term which is regarded as a contribution from orbit of zero-length, and it has peak at 
s = 0. The remaining part has a functional form exhibiting successive peaks at the periods 
of classical periodic orbits and their heights are proportional to the amplitude factors of 
the corresponding orbits. By comparing the calculated F'^'i™)(s) with F'^'"^^(s), we can thus 
extract information about periodic orbits from the quantum spectrum. 



4. SHELL STRUCTURE ENERGY CALCULATION 

A useful quantitative measure of shell structure is the shell structure energy which is 
defined as the fluctuation part of the sum of single-particle energies, i.e., 

N 

£,i,{N)=J2Ek-£{N), (4.1) 
fe=l 
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where E^'s represent eigenvalues of the single-particle Hamiltonian. The second term £{N) 
is obtained by smoothing the first term by means of the Strutinsky method; it is a smooth 
function of particle number N and of other potential parameters. The shell structure energy 
takes a large negative value when the single-particle level density at the Fermi surface is 
low. 

We have carried out a systematic calculation of the shell structure energy for the Hamil- 
tonian (3.1) as a function of the deformation parameters Jose (quadrupole deformation), A30 
(octupole deformation) and of the particle number N. Figure |l| shows a map of the lo- 
cal minima of shell structure energies in the two-dimensional deformation parameter space 
calculated for particle numbers in the range 16 < < 160. The centers of discs show 
the loci of local minima (corresponding to different values of N) and the sizes of the discs 
represent the absolute values of shell structure energies. As the order of magnitude of the 
shell structure energy is roughly proportional to N^^"^ in the harmonic oscillator case [ p^ , 
we normalize them by multiplying N^^^^ in order to compare system with different values 
of iV. 

For A30 — 0, it is well known that strong shell structure exists at Jose — (spherical), 
3/5 (prolate superdeformed) , 6/7 (prolate hyperdeformed) , —3/4 (oblate superdeformed) 
and so on. The distribution of discs is actually dense around such points. We find in Fig. ^ 
that, in addition to such known cases, prominent shell structure emerges also at finite values 
of octupole deformation. Their shell-structure energies are comparable, in magnitude, to 
(sometimes larger than) those for the purely quadrupole shapes. The most remarkable region 
is that of A30 = 0.3 ~ 0.4 and Sqsc — 0.1. Using the semiclassical theory, let us analyze in 
the following sections the mechanism which creates the above new shell structures for the 
combination of the quadrupole and octupole deformations. 



5. SEMICLASSICAL ANALYSES 

In this section, we investigate how the structure of quantum spectrum changes as the 
octupole deformation parameter A30 is increased, fixing quadrupole deformation parameter 
(5osc at several positive values (between the spherical and the prolate superdeformed shapes). 
Figure ^ shows the single-particle spectra as functions of A30 for the Hamiltonian ( |3.l[ ) with 
Sosc = 0.1, 0.3 and 0.5, The corresponding axis ratios u}±/uJz are 31/28, 11/8 and 7/4, 
respectively. They are not in simple ratios so that there is no prominent shell structure 
at A30 = 0. However, new shell structures emerge at finite values of A30. To see features 
of these spectra, we show in Fig. |^ the Fourier transforms of the level density defined by 



(3.10). As discussed in section 3 3.2, one see prominent Fourier peaks at the periods of 
classical periodic orbits. It means that the fluctuation of the spectrum is characterized by 
the periodic orbits, demonstrating a beautiful applicability to our model of the semiclassical 



method in section 2 2.2. For elucidating the features of shell structure, it is essential to 
understand the behavior of these Fourier peaks with respect to the deformation parameters. 
As discussed in the previous sections, short periodic orbits play important roles for the 
formation of gross structure of quantum spectrum. There are various periodic orbits of 
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various topologies in each part of the deformation parameter space and the same type of 
orbits change their characters as the parameters change. This fact is clearly seen in the 
Fourier transforms where peaks corresponding to certain orbits change their heights. The 
heights of the peaks represent nothing but the strengths of shell structures. Let us discuss 
in the foUowings what kind of periodic orbits exist and how they determine the features of 
quantum spectra in several regions on the ((5osc, A30) plane. 

Case of Sqbc = 0.1 

Let us first take up the case of 5osc = 0.1, where we obtain especially strong shell structures 
at finite A30 values in the shell structure energy calculation (see Fig. and let us discuss 
which orbits are responsible for these shell structures. Figure ^ shows, for several values 
of A30, some planar periodic orbits in the plane including the symmetry axis. We use the 
Monodromy Method developed by Baranger et al. ||l^ to calculate periodic orbits and their 
monodromy matrices. At A30 = 0, the most important orbit family is the ellipse-shaped one 
in the (x, y) plane. The next orbit family is 31:28 Lissageous. They are very long orbits and 
unimportant for the gross shell structure. Adding the octupole deformation, new types of 
orbits are born. Orbit PR appears at A30 = 0.12 by the isochronous bifurcation of the orbit 
PA, and PA becomes unstable after this bifurcation. At A30 = 0.24, the orbit PA becomes 
stable again and a new orbit PM appears. 

To see these bifurcations, the Poincare map is a very convenient implement. As our 
Hamiltonian has axial symmetry, we can treat it as a two-dimensional system with fixed 
angular momentum p^. Figure^ is the Poincare map {z,pz) for the 'projected' Hamiltonian 
( |3.5[ ) with Pip = 0. Here the surface is defined by pp = and pp < 0. In the top panel, 
the center of the main torus corresponds to the orbit PA. In the middle panel, the orbit PA 
becomes an unstable saddle and a new pair of island is created. These islands correspond 
to the orbit PR. In the bottom panel the orbit PA comes back to a stable one and creates 
a new pair of saddles confronting horizontally. These saddles correspond to the orbit PM. 

As is well known, stabilities of periodic orbits are determined by the monodromy matrices. 
Figure |^ shows the traces of the two-dimensional monodromy matrices for orbits associated 
with the above bifurcations. The bifurcation occurs at TrM = 2 where the monodromy 
matrix has the unit eigenvalues (1,1). Period n-upling bifurcations occur at TrM" = 2. 
Higher order (large n) resonances occur at every points of the torus and new orbits bifurcate 
from them, but they are of rather long periods and are related with more detailed structure 
of the spectrum. In the Fourier transform Fig. ^(a) , the peak corresponding to the orbit PA 
with period T ~ 2tt/uj± strongly enhances at A30 — 0.2. It mainly characterizes the shell 
structure seen in the spectrum Fig. ^(a) at A30 = 0.2 0.3. This enhancement is related 
with the above bifurcations of PA at A30 = 0.128 and 0.227. 
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Case of Sqsc = 0-3 

Next let us discuss the case of (5osc = 0.3. Some important classical orbits are drawn 
in Fig. 1^. Orbit PO is born out of the period-doubling bifurcation of the orbit IL at 
A30 = 0.215. Orbits PP and PQ are created by the saddle-node bifurcation (pair creation 
of stable and unstable orbits from nothing) at A30 = 0.221. Orbit PR emerges from the 
isochronous bifurcation of the orbit PA at A30 = 0.338. In the Fourier transform Fig. ^(b), 
the enhancement of the peak with s ~ 1 is related with the isochronous bifurcation. The 
orbits associated with the period-doubling and saddle-node bifurcations mentioned above are 
of similar periods and contribute to the same peak with s ~ 2.75 in the Fourier transform. 
It also shows strong enhancement in the bifurcation region. 

Case of (?osc = 0.5 

Lastly, let us consider the case of 6osc = 0.5. This quadrupole deformation is almost 
equivalent to the case of axis ratio \/3 : 1 treated in Ref. [|l8| . The feature of shell is rather 
weak at A30 = 0. However, one may notice in Fig. |^ that a remarkable shell structure 
emerges at A30 ~ 0.3. In the Fourier transform Fig. ||(c), a significant enhancement of the 
peak with s ~ 1.75 is observed. It seems that this is the most typical example in our model, 
which clearly exhibits a new shell structure emerging at finite octupole deformation. 

Some short periodic orbits are illustrated in Fig. ||. At A30 = 0, the ellipse-shaped family 
in the (x, y) plane characterize the structure of the spectrum. Orbit PB is born out of 
the isochronous bifurcation of the orbit IL at A30 = 0.276, orbits PC and PD are created 
by the saddle-node bifurcation at A30 = 0.274. The Poincare maps and the traces of the 
monodromy matrices are displayed in Figs. ^ and respectively. The strong enhancement 
of the peak with s ~ 1.85 at A30 — 0.4 are related with the bifurcations mentioned above. 
Orbits associated with these bifurcations are of similar periods and thus contribute to the 
same peak in the Fourier transform Fig ^(c). These bifurcations may be regarded as the 
mechanism which creates the prominent shell structure at finite octupole deformation. 

Orbits PE and PF are born associated with the period-tripling bifurcation of the orbit 
PA at A30 = 0.376. A quantum signature of this bifurcation is also seen in the Fourier peak 
with s ~ 3. Thus, we can conclude that the enhancement of Fourier peak in the bifurcation 
region is a general phenomenon. 

6. BIFURCATION LINES IN THE TWO-DIMENSIONAL PARAMETER SPACE 

In the preceding section, we presented several examples where bifurcations of short pe- 
riodic orbits play important roles in forming shell structures. In this section, let us dis- 
cuss roles of the bifurcations varying two deformation parameters systematically. In the 
n-dimensional parameter space, bifurcation points generally form (n— l)-dimensional man- 
ifolds; they are 1-dimensional curves in the present case. We have evaluated these bifurca- 
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tion lines for some short periodic orbits. The result of calculation is presented in Figure |T^. 
These bifurcation lines always emerge at the points where the ratios uj^/lOz are rational 
and A30 = 0. Note that (5osc=0, 3/5, 6/7 and —3/4 correspond to the spherical, prolate 
superdeformed, prolate hyperdeformed and oblate superdeformcd shapes, respectively. The 
orbit PA causes isochronous bifurcations along the lines i and ii (in direction from the 
right to the left). As discussed in the case Jose = 0.1, new orbits PM and PR are born 
after the bifurcations i and ii, respectively. On the other hand, as discussed in the case 
5osc = 0.5, the lines iii'-^^vi are related with the orbits PB, PC, PD, IL and 2PA (double 
traversals of the orbit PA). The orbit IL causes isochronous bifurcation (from the left to 
the right) and produces the orbit PB along the line iii, while the orbit PC (unstable) 
and PD (stable) are born out of the saddle-node bifurcation along the line iv in the region 
Sqsc ^ 0.5. As SoBc increases, orbits PB and PC exchange their stabilities with each other 
along the line v and annihilate into 2PA along the line vi. Along the line vii, the orbit PA 
causes period-tripling bifurcation (from the right to the left) and produces orbits PE and 
PP. 

According to the discussions in the preceding sections, we expect strong shell effects along 
these bifurcation lines. Prominent shell structures are known to exist for the spherical and 
superdeformed shapes (the end points of the bifurcation lines) where periodic orbit condi- 
tions for the harmonic oscillator potential are satisfied. Conditions for the emergence of 
pronounced shell structure in non-integrable systems are not well known, however. As is 
indicated in Fig. |l|, shell structure energies at equilibrium shapes with finite A30 are compa- 
rable in magnitude to those of spherical and superdeformed shapes. The most remarkable 
one is the region with (5osc — 0.1 and A30 = 0.3 ^ 0.4. The bifurcation map of Fig. |ll| 
suggests that the strong shell effect in this region may be connected with the bifurcation 
of orbit PA (the lines i and ii). In fact, we saw in the preceding section that the Fourier 
peak corresponding to the orbit PA is strongly enhanced by the bifurcation phenomena. It 
is worth emphasizing that such a strong shell effect can arise associated with classical orbit 
bifurcations in non-integrable systems. 

One should also note that most of the large discs in Fig. |l| at finite A30 locate in the 
prolate side (Jose > 0). This is related with the property of the shortest periodic orbit: 
For prolate shapes, the shortest orbit is of the type-PA. It is degenerate and stable, so its 
contribution to the level density is important. On the other hand, the shortest orbit for 
oblate shapes is of the type-IL, which is isolated and, accordingly, its contribution to the 
level density is rather small. The orbit of type-PA, the second shortest orbit, is unstable 
against the octupole deformation, and less important in comparison with the prolate case. 
A more detailed discussion on the difference in stability between the oblate and prolate 
superdeformed shapes against octupole deformations will be given in the next section. The 
same problem was discussed also in [ pO|j2H | from somewhat different point of view. 
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7. OCTUPOLE DEFORMATION SUPERPOSED ON THE PROLATE AND 
OBLATE SUPERDEFORMATIONS 

Let us discuss the origin of the difference in octupole stabihty between the prolate and 
oblate superdeformed states. In Ref. fl^jlj] we have shown that the supershell effect in the 
prolate superdeformed states increases with increasing octupole deformation. As an under- 
lying mechanism of that enhancement, we emphasized the importance of stability properties 
of two kinds of periodic orbit family and of their interference effect. The oblate case is sim- 
ilar to the prolate case in that there are two kinds of periodic orbit family whose periods 
are in the ratio 2:1. But the structure of quantum spectrum is quite different with each 
other. In Fig. |l^ are compared the single-particle spectra for the prolate {uj±/u}z = 2) and 
oblate {u}±/ujz — 1/2) superdeformed oscillators as functions of the octupole deformation 
parameter A30. The way the degeneracy be solved is different between the two. The oc- 
tupole operator Y30 has matrix elements between states in the same major shell in the oblate 
case and, therefore, it affects the spectrum in the first order perturbation, while it affects 
only in the second order in the prolate case. Let us discuss below how this difference be 
explained in terms of the classical dynamics point of view. For this purpose, representative 
periodic orbits with short periods are displayed in Fig. ^ for several octupole deformation 
parameters A30. 

We first compare the features of the two spectra without the octupole term. Figure 
shows the oscillating level density smoothed to an energy width SE = uJsh/2 (wsh being 
the energy spacing between adjacent major shells). A characteristic feature is that the 
prolate superdeformed spectrum has an undulating pattern (supershell structure) while 
oblate one does not. This is due to the difference in degeneracies of contributing periodic 
orbits. As discussed in Ref. [Q, two orbit families (corresponding to orbits (PB, PC) 
and orbit PA in the upper panel of Fig. respectively) in the prolate superdeformed 
states have degeneracies 4 and 2. In the oblate case, the orbit family with period 2tt/uj±^ 
(corresponding to orbits (PA, PL) in the lower panel of Fig. ^ has the maximal degeneracy 
4, but the shortest orbit (the linear orbit IL along the z-axis) is isolated and has degeneracy 
1. Thus the interference effect between these two families is so small that one cannot see 
the supershell effect in the spectrum. 

Figure ^ compare the Fourier transforms of the level density for the prolate and oblate 
cases. One can hardly see the component at s = 1 in the oblate case and the oscillating pat- 
tern of the spectrum is determined by the s = 2 component almost exclusively. Comparing 
the two figures, one notice that the reduction rate of the peak-height due to the octupole 
deformation is much greater in the oblate case. This rapid decline clearly corresponds to 
the rapid disappearance of the shell effect in the oblate case. 

The main reason for reduction of the shell effect with increasing octupole deformation 
is two-fold: The first is the reduction of degeneracy of the periodic orbit famihes, and the 
second is the change of stabilities. As the degeneracies are the same for the major orbit 
families in both cases (orbits PB, PC, • • • in the prolate case and orbits PA, PL, • • • in 
the oblate case), we expect that the differences are associated mainly with the stability 
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properties. In Fig. ^ we show the stabihty factors X = ^J\det{l^^JAr)\ calculated as 
function of A30. The stability factors for orbits in the oblate potential depend linearly on 
A30 for A30 — while they depend quadratically in the prolate case. Consequently the 
amplitude factors reduce much faster in the former case. This seem to be the main cause 
of the rapid reduction of shells with octupole deformation in the oblate potential. 



8. SUMMARY AND CONCLUSION 



We have analyzed the gross structure of single-particle spectra in reflection-asymmetric 
deformed oscillator potentials using the semiclassical method. Our model is non-integrable 
and is regarded as a mixed system where regular and chaotic dynamics coexist. The periodic 
orbit theory, which is well established for regular and strongly chaotic limits, seems to be 
also applicable to such a situation. Fourier transforms of the quantum level density reveal 
almost perfect correspondences with classical periodic orbits. Importance of classical orbit 
bifurcations has been demonstrated in our model. Strong shell effects arise also for rather 
chaotic regions, and their strengths are comparable in magnitude to those of regular regions. 
We obtain an interesting result which indicates that classical bifurcations may be responsible 
for the emergence of shell structure in the mixed system. Applications of the semiclassical 
theory of shell structure to more realistic mean-field potential models and identifications 
of classical orbits which play decisive roles in determining exotic shapes of nuclei or micro- 
clusters remain as exciting future subjects. 
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vmin.ps 

FIG. 1. Local minima of shell structure energy in the two-dimensional deformation parameter 
space (i5osc,A3o). Size of each disc represents the absolute value of the shell structure energy 
normalized as Ssi^/N^^^ . Plotted are for even A'^ in the range 16 < A'^ < 160. 
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FIG. 2. Single-particle spectrum of the 



Hamiltonian (3.1) with deformation parameter 
Sosc = 0.1, 0.3 and 0.5 as functions of the oc- 
tupole parameter A30. Dashed and solid curves 
represent the levels whose K quantum numbers 
are zero and nonzero, and the latter degenerate 
in two due to the time-reversal symmetry. 



(a) i5osc = 0.1 
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(b) Sosc = 0.3 
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Fig. y, to continue. 


(c) (5osc = 0.5 
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FIG. 3. Fourier transforms |-F(s)| of the level densities for the Hamiltonian (3.1) with Sosc ~ 0.1, 
0.3 and 0.5, plotted as functions of the action s and the octupole parameter A30. 
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FIG. 4. Some short planar periodic orbits of the Hamiltonian (3J) with Sosc =0.1 
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FIG. 5. Poincare surface of section {z,pz) for the Hamiltonian (3.5) with = 0. The surface 
is defined by pp = and pp < 0. 
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FIG. 6. Traces of the Monodromy matrices at Sosc = 0.1 for some 
short periodic orbits plotted as functions of the octupole deformation 
parameter A30. 
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FIG. 7. Some short planar periodic orbits for the Hamiltonian (3.1) with Sosc = 0.3 
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FIG. 8. Some short planar periodic orbits for the Hamiltonian (3.1) with Sosc = 0.5 
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FIG. 9. Poincare surfaces of section for Sosc = 0.5 in the bifurcation region of short periodic 
orbits. The surface {p,Pp) is defined by 2: = and pz > 0. 
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FIG. 10. Same as Fig. H but Sosc = 0.5. 
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FIG. 11. Bifurcation lines for some short periodic orbits in the two-dimensional deformation 
parameter space (5osc,A3o). 
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(a) (b) 
FIG. 12. Single-particle spectra for (a) prolate and (b) oblate superdeformed oscillators plotted 

as functions of the octupole deformation parameter A30. 
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FIG. 13. Oscillating level density of prolate and oblate superdeformed os- 
cillators, which are smoothed to an energy width SE = hush/'i. 
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FIG. 14. Some short periodic orbits in the prolate (upper panel) and the oblate (lower panel) 
superdeformed potentials with octupole deformations. 
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FIG. 15. Fourier transforms of the quantum level densities for the prolate (upper panel) and 
the oblate (lower panel) superdeformed oscillators with octupole deformations. 
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FIG. 16. Stability factors X = det(l — M,-) of the short periodic orbits for the prolate (left) 
and the oblate (right) superdeformed potentials plotted as functions of the octupole deformation 
parameter A30. The sign of X is that of det(Mr — 1) = Tr A/,- — 2. Namely, —2 < X < for stable 
orbits and otherwise for unstable ones. 



